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ABSTRACT 

The Carnegie Hubble Program (CHP) is designed to improve the extragalactic dis- 
tance scale using data from the post-cryogenic era of Spitzer. The ultimate goal is a 
determination of the Hubble constant to an accuracy of 2%. This paper is the first 
in a series on the Cepheid population of the Large Magellanic Cloud, and focusses on 
the period-luminosity relations (Leavitt laws) that will be used, in conjunction with 
observations of Milky Way Cepheids, to set the slope and zero-point of the Cepheid 
distance scale in the mid-infrared. To this end, we have obtained uniformly-sampled 
light curves for 85 LMC Cepheids, having periods between 6 and 140 days. Period- 
luminosity and period-color relations are presented in the 3.6 /im and 4.5 /im bands. 
We demonstrate that the 3.6 /jm band is a superb distance indicator. The cyclical 
variation of the [3.6]-[4.5] color has been measured for the first time. We attribute the 
amplitude and phase of the color curves to the dissociation and recombination of CO 
molecules in the Cepheid's atmosphere. The CO affects only the 4.5 /jm flux making it 
a potential metallicity indicator. 



Subject headings: Cepheids — distance scale — infrared: stars — Magellanic Clouds 



- 2 - 



Introduction 



In the past decade, t he "factor-of-twq cont r oversy" surrounding the value of the Hubble 



constant has been resolved (Freedman et al. 



Riess et al. 



2011 



2001 



Spergel et al 



2007 



Freedman &; Madore 



2010 



). However, the uncertainty on Ho currently remains at the level of a few percent 
and dominates the error budget when deriving other cosmological parameters such as the equa- 
tion of state parameter, ujq. Reducing the uncertainty on H n wi l l enab le us to further constrain 
the parameters in the concordance cosmological model — as iHul (|2005l ) states: "The single most 
important complement to the CMB for measuring the dark energy equation of state at z ~ 0.5 is 
a determination of the Hubble constant to better than a few percent." 

The Carnegie Hubble Program (CHP) is a "Warm Spitzer" legacy mission (Spitzer Exploration 
Program PID 60010 "The Hubble Constant", P.I. Freedman). It has the primary goal of reducing 
the systematic uncertainties on the Hubble constant to 3% using Spitzer data alone, an d ultimately 



2% w ith the addition of JWST. The program uses the 3.6 and 4.5 fim bands of IRAC (jFazio et al 



20041 ) to recalibrate the Cepheid distance scale in the mid-infrared, moving out in to the Hubble 
flow by calibrating and applying a revised mid-IR Tully-Fisher relation. 



A detailed description of the Carnegie Hubble Program is given in iFreedman et al.l (|201ll ) ; an 
overview of our strategy follows below. 

The zero-point of the Cepheid period-luminosity (PL) relation (Leavitt law) will be calibrated 
using Galactic Cepheids. This sample contains the 10 Cepheids with tri gonometric parallaxes mea- 
sured using the Fine Guidance Sensors on the Hubble Space Telescope ([Benedict et al.l 120071 ). We 
also observed 27 Cepheids within 4 kpc of the Sun — close enough for GAIA parallax measure- 
ments. Seventeen of these are known members of Galactic open clusters for which main sequence 
fitting distances are available. The Milky Way program is described in Monson et al. (2011, in 
preparation) . 

The second rung on the distance ladder (the subject of this work) involves the long-period 
(P > 10 days) end of the Leavitt law. We will calibrate the slope and dispersion of this period 
range using a sample from the Large Magellanic Cloud. The LMC is an important stepping stone 
for the extragalactic d istance scale; it was used as the base for the Key Project distance scale 



( Freedman et al 



2001 



but uncertainty in its distance eventually became the dominating system- 
atic. The calibration found here will then be applied to nearby galaxies in the Local Group, forming 
a Cepheid distance ladder based on a single telescope and instrument. Additionally, by including 
galaxies such as M33, the effect of metallicity on the PL slope and zero-point will be studied (e.g., 



Scowcroft et al.l 12003 ). The calibration will then be applied at a step further out, extending to 



galaxies beyond the Local Group, such as M81 and NGC 4258, the well known maser galaxy. 

Next, we will extend our distance ladder into the Hubble Flow using the Tully-Fisher relation, 
itself re-calibrated in the mid-infrared. The dispersion of this relation is expected to drop at 
IRAC wavelengths, simply due to the decreased sensitivity to inclination-dependent extinction 



- 3 - 



corrections, and will be calibrated using galaxies that contain Cepheids. The distances to over 500 
target galaxies in our program will be measured. Finally, we will complete our distance ladder 
by adding 44 supernova-host field galaxies. These steps will allow us to measure accurate and 
consistent distances from the closest Cepheids to the most distant galaxies. 

This work deals with Cepheids in the Large Magellanic Cloud (LMC). The LMC Cephe i d pop - 
ulation has been studied extensively for over half a century. For example, IShaplev &: Naill ()1940l ) 
presented well-sampled light curves for 40 Cepheids in the L MC, and study of the p opulation con- 
tinues to this day with large-scale projects such as OGLE III ( Soszvnski et al.ll2008l ). Its proximity 
means we can readily observe Cepheids over a range of periods, from those that are similar to many 
of the Galactic parallax sample (i.e. P < 10 days), to the long period Cepheids with P > 10 days, 
which more generally overlap Cepheids most easily observed in our most distant targets. 

Until relatively recently, the majority of distance measurements have been undertaken at op- 
tical wavelengths (e.g. F01). Cepheid studies at optical wavelengths have their shortcomings, 
the main one being extinction. Although the effect can, to first order, be removed by use of the 
reddening-free Wesenheit index (|Madordll976l ) this technique still requires prior knowledge of the 
extinction law, as well as an assumption that it is universal. By mo ving to the mid-infrared , 
reddening and extinction are diminished by around a factor of twenty (jRieke fc Lebofskylll985l ). 
making their absolu te contribution and thei r uncertainties negligible. In addition to the drop 
in extinction effects, iMcGonegal et al.1 ()1982l ) describe two other advantages of the infrared over 
the optical. In the infrared the amplitudes of the Cepheids' lightcurves decrease, as does the in- 
trinsic width of the instability str i p, be cause these wavelengths are less sensitive to temperature 
changes. Indeed IMcGonegal et al.l (|1982l ) demonstrated that the width of the near-infrared H band 
period-luminosity relation from random observations is less than the width from time-averaged B— 
band observat i ons. Near-infrared observations of LMC Cepheids were more recently obtained by 
Persson et al.l (12004! ) . By moving to longer wavelengths in the mid-infrared, in combination with 
well phased observations, we can decrease the measured width even further. 

The IRAC imager on Spitzer is a superb instrument for undertaking a recalibration of the 
Cepheid distance scale. Like Hubble, Spitzer has the advantage of operating in a stable environment 
without weather or seeing variations, and with great flexibility in scheduling. As such, we have 
been able to obtain precise and deterministically well-sampled light curves for 85 Cepheids. The 
observations of the LMC Cepheids were one of the first programs undertaken by post-cryogenic 
"Warm Spitzer". 

The target selection and observations are described in Section [2] and the photometry and 
calibration are discussed in Section [3j Light and color curves are presented for each Cepheid in 
Section^ PL relations, including a discussion on their use in determining the tilt of the LMC, are 
given in Section [5j The period-color relation at mean light has been measured for the first time at 
these wavelengths; it is discussed in Section [6l Section [7] provides a summary. 
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Observations 



2.1. Target Selection 



The LMC has a large sample of well-studied C epheids to choose from, as light curves have been 
measured for many thousan ds. The OGLE studies (jUdalski et al, l ll992l : lUdalski. Kubiak fc Szvmanskil 
19971 ; ISoszvnski et al.ll2008l ). while having a primary goal of observing micro- lensing events, have 
produced extremely large Cepheid catalogs for th e LMC and SMC in t he optical. For this study, we 
have chosen a subset of the sample observed by iPersson et al.l (|2004l . hereafter P04), who studied 
92 LMC Cepheids in the infrared JHK S bands. These Cepheids also had optical (UBVRI) light 
curves available. The majority of them are free from crowding, making them ideal for this study. 
We chose the 85 Cepheids with periods greater than about 6 days. This ensures that the period 
distribution overlaps with both the Galactic sample, as well as the most remote galaxies in our 
program. The Cepheids are distributed over the face of the LMC, allowing any known inclination 
effects (c.f., P04) to be re-measured and accounted for. 



2.2. "Warm Spitzer" Observations 

The data were obtained in programs 61000, 61004, 61005, 61006, and 61007, executed between 
2009 October 3 and 2010 July 18, at the start of the "Warm Spitzer" mission. Each Cepheid was 
observed at 24 epochs. For those Cepheids with P > 12 d, the observations were taken over a 
single cycle, spaced on average by P/24 days. This produced regularly-sampled light curves while 
preserving the high scheduling efficiency of the spacecraft. The shortest-period Cepheids (P < 
12 d) could not be observed in such a tightly-phased manner; instead we scheduled 24 observations 
separated by 12 ± 4 days. The Appendix explains the choice of 24 evenly-spaced epochs. 

Each IRAC observation consisted of ten frames, five in each of the 3.6 (im and 4.5 /im bands, 
with a medium-scale, 5-point gaussian dither pattern. (Dithering of exposures is essential for 
accurate photometry.) The frame time was 2.0 s, giving an effective per-pixel integration time of 
1.2 s in each frame. IRAC is known to be susceptible to persistent images (i.e. latent images of 
bright, but often unsaturated, stars can remain in later frames) and this effect is apparently worse 
in Warm Spitzer than in the original mission (see IRAC handbook, version 2.Q£|, section 7.2.8). 
By using a large scale dither pattern, where subsequent observations of the Cepheid are spaced by 
approximately half the chip, we have mitigated the possibility of our photometry being affected by 
previous observations. 



: http: //irsa. ipac . caltech.edu/data/SPITZER/docs/irac/ iracinstrumenthandbook/ 
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3. Data reduction and Photometry 

These early "Warm Spitzer" data were some of the first to be reduced with the new "Warm 
IRAC" pipeline. The images were processed using the Spitzer Science Center pipeline version 
S18.18.0, which does not include the artifact correction steps present in the cold mission pipelines. 
Therefore, in addition to the normal Spitzer Science Center (SSC) processing, the BCD (Basic Cal- 
i brated Data) frames we re run through the Warm-Mission Column Pulldown CorrectoJE mopex 



(jMakovoz fe Khanll2005l ) was used with the standard overlap and mosaic name lists for the overlap 



correction and mosaicking. 

The most important objective of this work is to produce an accurate calibration of the Cepheid 
period-lu minosity rela t ion. Our photometry is on the same scale as the standard calibration de- 



scribed in lReach et al.l (|2005l ). They define the flux scale and zero point such that they correspond 
to a 10-pixel radius aperture, with a sky annulus running from 12 to 20 pixels (in units of native 
IRAC pixels) with pixel-phase and array-location corrections applied as necessary. 

Although the majority of the Cepheids in our sample are isolated, several have close neighbors 
that fall inside the 10-pixel aperture. These stars could not be measured correctly using aper- 
ture photometry. In order to reduce the entire LMC sample in a consistent way, Point Response 
Function (PRF) fitting was used in place of aperture photometry, utilizing the SSC apex software 



(jMakovoz Sz Marleaull2005l ). The point response function (PRF) models the detector response to a 
point source, and can be thought of as an over-sampled representation of the point spread function. 
The BCDs were mosaicked to obtain the bad pixel maps, but the PRF fitting was performed on 
each BCD individually using the apex-lframe name list, with two changes. First, rather than using 
the PRF corresponding to the center of the array, the set of 25 spatially-dependent PRFs was 
used for each channel. The IRAC detectors are known to have some spatial response variations; 
using the PRF map ensured that any variation in flux was properly accounted for. Second, the 
normalization radius - the radius to which the fit is calculated - was set to 1000 pixels. As the 
PRF is sampled at 100 times per pixel, the 1000 pixel normalization radius corresponds to the 
calibration aperture radius. For each BCD, all the stars in each image were measured. 



3.1. PRF calibration 

The construction of the IRAC PRF is a complex task requiring a large number of observa- 
tions of a stable star. Warm Spitzer PRF models are not yet available from the Spitzer Science 
Center; however, it is still possible to obtain accurate photometry using the cold PRF, correcting 
as necessary for the predicted differences between the cold and warm models. The most significant 



http: //ssc . spitzer . caltech. edu/dataanalysistools/tools/ contributed/irac/f ixpulldown/ 



3 PRFs are available from http://ssc.spitzer.caltech.edu/irac/calibrationfiles/psfprf/ 
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difference between the two models is the pixel-phase correction (PPC). 

The PPC is necessary given that the response of the IRAC detectors is known to vary on a 
sub-pixel level. This means that the measured flux of an object not only depends on the position 
of the star in the frame, but also on the position of the peak of the PSF within the pixel. The 
fractional position within a pixel is known as the "pixel phase". In the cryogenic mission the 
pixel-phase effect was around 4% (peak-to-peak), but is approximately double that in the warm 
mission. In addition to the change in magnitude, whereas the cold mission PPC was represented by 
a symmetric Gaussian centered on the peak response position, the warm mission PPC is described 
by a Gaussian, differing in width in the x and y directions. 

As the PRF is an oversampled representation of the pixel's response to light, it includes the 
cold mission PPC by definition. To correct the warm mission data, the cold PPC must first be 
multiplied out of the apex photometry results, and the warm correction then divided in. Although 
it is possible that the pixel-phase effect will cancel out when the results from a dither pattern are 
averaged together, in our data the PPC was applied to each BCD before the averaging took place. 



3.2. Photometric stability 

The stability of IRAC is of fundamental importance to our program. We have quantified the 
stability of the instrumental zero-point by measuring constant stars in our Cepheid fields, and 
monitoring their magnitudes over time. The photometry method was identical to that used on 
the Cepheid targets. In Figure [1] we show the photometry for constant stars in two fields. Star 
A was chosen in the HV00872 field, at a = 4 /l 55 m 17.09 s , 6 = -67°27'53.82" , and Star B was 
in the HV02883 field at a = 4 h 56 m 21.24 s , 5 = -64°41'35.65" . In both cases, the dispersion of 
the measurements is approximately ±0.01 mag, and no systematic trend with time is seen. The 
Cepheids in the fields have very different periods (29.82 d vs. 108.97 d), hence the time baseline 
and the cadence of the observations differ. This confirms both the short and long term stability of 
IRAC. 



3.3. Photometry comparison 

As mentioned above, the IRAC magnitude system is defined by the flux measured within a 
10-pixel aperture. To verify that the photometry for this work is on the calibrated system, the data 
for the Cepheid HV00872 were also processed using the SSC IDL routine photo Jorstar.pro (SSC 
2010, private communication). This routine performs aperture photometry on a single star using 
two aperture sizes — a 3 pixel aperture with a 3-7 pixel sky annulus, an d the 10 pixe l aper ture 
with a 12-20 pixel sky annulus used in the calibration definition given by iReach et al.l (|2005l ). It 
also includes the latest pixel-phase corrections and array-location corrections. 
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The APEX photometry was compared to the results from both aperture sizes. No systematic 
trends with magnitude were found. The difference between the apex and the 3-pixel aperture 
fluxes were consistent with the published aperture correction for transformations between the 3 
and 10 pixel apertures. Both the apex and small aperture results were then compared with the 
large aperture values. As was expected, the scatter in the large aperture results was much larger 
than the other two methods, but no systematics were found. Aperture corrections were derived 
for each band, and agreed with the cold mission values (1.124 and 1.127 at 3.6 /xm and 4.5 /im 
respectively) to within 1%. 

To put the PRF photometry on the same system as the standard aperture photometry the 
following correction was used: 

_ F PRF xAx PPCpoid ^ 
PPC W 

arm 

where F is the flux in the standard system, Fppp is the flux output by apex, A is the aperture 
correction (1.03 and 1.02 for channels 1 and 2 respectively) and PPCcoid an d PPCwarm are the 
cold and warm pixel phase corrections. The fluxes were converted to magnitudes using the standard 
zero-magnitude flux densities of 280.9 Jy ([3.6]) and 179.7 Jy ([4.5]). 



Results 



4.1. Periods 



We first refined the available peri ods as follows. Light curves were produced in the 3.6 and 
4.5 £tm bands. The data points from iMadore et al.1 (|2009l ). taken as part of the SAGE project 
(jMeixner et al.ll2006l ). were also phased in. The periods given in P04 were taken as initial values; 
these were improved by iterating on the period and re-phasing the data until a minimum dispersion 
in the light curves was obtained in all available wavebands (ET, B,V, R, I, J, H, K, [3.6] , [4.5] , [5.8] , [8.0] , 
but note that not all wavelengths were available for all Cepheids). On average, the periods changed 
by less than 1%. In some cases it was clear that the period of the Cepheid had changed in the 
time interval between the optical and infrared observations. This is to be expected, as we utilized 
as much archival data as possible; the time baseline for some Cepheids was of the order of many 
decades, which is comparable to the time many Cepheids have been found to signifi cantly change 
their periods (for a discussion of period changes in a sample of LMC Cepheids, see iPietrukowicz 



2001 



). The practice of fitting all wavebands simultaneously makes it easy to find objects for which 
this is the CclSG, clS the period that fits the later observations will result in an incorrectly phased 
light curve for the archival data, and vice versa. When this was the case the IRAC and JHK data 
alone were used to refine the period and to phase the data. The finally adopted periods are given 
in Table [TJ 
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4.2. Infrared Light Curves and Mean Magnitudes 

Once the periods had been adjusted, the IRAC light curves were plotted using just the CHP 
data. This avoids issues such as slight period changes and larger photometric uncertainties in the 
archival data that could affect the IRAC fit. The resulting light curves are shown in Figure [2] in 
order of increasing period. The GLOESS program (described in section 3.2 of P04) was used to 
interpolate between the data points and to thus draw the smooth curves. GLOESS was also used to 
determine time-averaged (intensity-mean) magnitudes and uncertainties. The results are given in 
Table [E 

For the [3.6] (top) and [4.5] (middle) panels, the statistical uncertainties on the individual 
measurements are comparable to the size of the plotted points. Table [2] gives a sample of the 
individual data points for HV00872; all data for this program are available with the online version 
of the article. 

The bottom panels of the plots in Figure [2] show the [3.6] — [4.5] color variation through the 
Cepheid pulsation cycle. As light curves of this quality have not previously been available in 
the near-IR, this cyclical effect has not been noted before. The effect is almost c ertainly due to 



absorp tion by the 4.65 fim CO bandhead in the 4.5 /jm filter, as first suggested bv lMarengo et al 



(2010a|). The color variations are discussed below in Section 



5. Mid- infrared Period Luminosity Relations 

Period-luminosity (PL) rel ations for the two IRAC ban ds were fitted using the weighted-least- 
squares technique described in lAkritas &: Bershadyl (119961 ). This takes into account the intrinsic 
width of the PL relation, giving a more robust result than a regular weighted-least-squares fit. 

The PL relations take the form 



M = o(IogP- l) + 6 



(2) 



and were fit using four different period cuts. The results of the different fits are given in Tab le [3] 
The relations derived from the two phase points available from SAGE (jMadore et al.l [2009J) are 
also given for co mparison. In Table [3 j N is the number of Cepheids in the fit, a and b were both 
found using the lAkritas Bershadvl (| 19961 ) weighted-least-squares method, a a and o\, are their 
respective uncertainties, and a is the standard deviation of the fit. 

The finally adopted relations are plotted in Figure O The solid line shows the adopted PL 
relations; these were the fits that used only Cepheids in the range 10 d < P < 60 d. Note the 
systematic deviation of the longest-period Cepheids, which are well known to be significantly fainter 
than the PL relation at virtually all wavelengths. Consistent with other studies (such as F01) we 
have excluded objects with P > 60 days from the fitting. The SAGE PLs (given in Table [3J are 
entirely consistent with the values found in this work. Although we find no evidence in our data to 
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suggest a break or change in slope in the P = 10 days region (c.f., ISandage. Tammann &; Reindl 



( 2004 )). we have also cut the low end of the period distribution, again consistent with F01. 



In both cases, the dispersion around the relation has fallen dramatically in comparison with 
previously published determinations. The previous measurements from SAGE gave standard devi- 
ations of 0.135 and 0.141 mag in [3.6] and [4.5] respectively, whereas the sca tter around our fina l 



relations are 0.108 and 0.115 mag. This is because the magnitudes used in the lMadore et al.l ((2009) 
PL relation are not fully phase-averaged, but rather they are the simple mean of two random phase 
points (see the Appendix). 

The precision of our phase-averaged light curves is best demonstrated by comparing the resid- 
uals from the [3.6] PL relation to those from the [4.5] relation. This is shown in Figure HI In this 
figure the small inset plot shows the entire data set, with the main plot showing just the objects 
with A mag < ±0.2, which encompassed 80 of the 85 Cepheids in the sample. It is clear that the 
residuals are highly correlated. That is, if a star lies above the mean PL relation in one band, it 
will tend to be above the mean in the other band as well. The same effect in the J and K$ bands 
was found by P04. 

Excluding two outliers, the range in A mag is ±0.2 mag. This spread is a combination of two 
effects: depth and tilt effects in the structure of the LMC, and the intrinsic magnitude width of 
the instability strip. The residual dispersion in this correlation is only 0.02 mag, suggesting that 
we are seeing the photometric limit of the data. 

Figure [5] shows how the residuals of the ([3.6] — [4.5]) color correlate with the PL residuals 
of the [3.6] magnitude. As we know from the previous plot that any dispersion in Figure [5] arises 
from the intrinsic properties of the Cepheids rather than photometric uncertainties, the width of the 
instability strip can be measured from this figure. Assuming that the instability strip is represented 
by a rectangular distribution, i.e. it is uniformly filled and has hard limits at the blue and red edges, 
the dispersion around the mean is related to the width by the relation a = R/\/V2, where R is the 
full width of the distribution. The color dispersion of the data is ±0.018 mag, leading to a color 
width of the instability strip of 0.062 mag. 



5.1. Tilt of the LMC 

By combining our knowledge of the Cepheid's location within the LMC with its position in 
the period-luminosity plane, the tilt of the galaxy with respect to the plane of the sky can be 
constructed, as was done in P04. 

In Figure [6] we show the positions of the Cepheids within the LMC in projection on the sky. 
The bottom left panel shows the position of each of the Cepheids relative to the center of the LMC. 
Solving for PL coefficients a and b, but now including £ and rj (distances from the center of the 
LMC in RA and Dec) as free parameters, we can obtain the tilt and the angle of the line of nodes if 
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the Cepheids on average lie in a plane. The top left and bottom right plots show the distribution of 
the points when viewed from orthogonal directions in the plane of the sky. Empty circles represent 
the Cepheids that are brighter than the average (i.e. they lie above the mean PL relation), and 
filled circles show those that are fainter than the average (below the mean PL relation). In this 
figure the reference coordinate system has been rotated by 50° so that the line of nodes of the LMC 
is parallel with the x-axis. In this plot, it is clear that the Cepheids in the bottom half tend to 
be fainter than the PL, hence further away, while those in the top half tend to be systematically 
brighter, hence closer than average. The simplest interpretation is that the LMC is in fact tilted 
with respect to the plane of the sky. 

By tilting the LMC by 28° relative to the plane of the sky we can minimize the width of the 
histogram in the top right corner of Fig. El This is in good agreement with the results of P04, and 
will be the subject of a more complete discussion of the data at a later time. 



6. The Role of CO Absorption 



We believe that the variations in the [3.6] — [4.5] color, both at mean light and throughout 
the pulsation cycle, are due to CO absorption in the 4.5 /iin band. We do not see any evidence of 
extended emiss i on sur roundi ng the Cepheid s in ou r images. This is consistent with the results of 
Marengo et aD (|2010al lbh and [Barmbv et al.l (|201ll ) who found evide nce for circumstellar em ission 
in the longer wavelength IRAC bands, but nothing at 3.6 or 4.5 ^m. iMarengo et alj (|2010al ) notes 
that the [3.6] — [4.5] color of the Cepheids is indistinguishable from their non-variable control stars, 
and is therefore caused by an intrinsic feature of stars of this spectral type, regardless of whether 
or not they are Cepheids. This rules out circumstellar envelopes and mass-loss as the source of the 
feature, leaving only the CO bandhead. 



6.1. Period Color Relation 

The possibility of a period- color (PC) relationship at mid-infrared wavelengths was correctly 



noted by IMarengo et al.1 (|2010al ). However, their magnitudes (and therefore colors) were measured 
at only a single epoch, introducing large scatter into the PC relation, as is shown in their figure 
5. Figure [7] shows the relationship between period and [3.6] — [4.5] color for our sample. The 
period-color relation was fit using Cepheids with 1 < logP < 1.8, and is found to be 

[3.6] - [4.5] = -0.087(±0.012)(log P - 1) + 0.005(±0.005). (3) 

The LMC Cepheids clearly exhibit a period-color relation in this period range. The standard 
deviation around the relation is ±0.018 mag. The relation appears to invert for the longest — 
period (logP > 2) Cepheids. Ignoring the logP > 2 sample, the long period Cepheids are brighter 
and bluer in [3.6] — [4.5] than their short period counterparts. At longer periods the Cepheids are 



-11- 



at cooler average temperatures, and so the carbon and oxygen are situated in CO molecules. As 
discussed previously, the CO absorbs some of the 4.5 fim flux, making the star appear bluer. As 
we move to shorter periods the average temperature rises, the CO is gradually dissociated, and 
the average colors become redder. The period-color relation reverses at P > 100 days. We do not 
know what causes this reversal, but the Cepheids at these periods are deviant in their luminosities 
and colors at almost all wavelengths. 



6.2. [3.6] - [4.5] Color Curves 

The period-color relation is closely related to the variations in color shown in Figure[2l as they 
are both due to CO affecting the 4.5 ^m flux. In a static stellar atmosphere, with a temperature 
and surface gravity corresponding to an F-type supergiant, the 4.5 fim flux is suppressed. In the 
case of Cepheids, however, we are not dealing with static atmospheres. As the Cepheid contracts 
its temperature increases. As the star becomes hotter, CO in the atmosphere dissociates, causing 
less absorption at 4.5 /um, making the star then appear redder than average. Once the Cepheid 
begins to expand again its temperature drops, allowing the C and O to recombine. As more CO 
is present, the 4.5 fim flux is once again suppressed, making the star bluer, in the particular color 
combination ([3.6] — [4.5]) discussed here. Inspection of Figure [2] shows that the phasing of the light 
and color curves exhibits precisely this behavior. 

There also appears to be a relation between period and the amplitude of the color variation, 
such that the short period Cepheids (P < 10 days) tend to have flat color curves. The cyclical 
variation appears to "turn on" around P = 10 days, increasing in amplitude with period. Again, 
this is a temperature effect. The short period Cepheids are so hot that the CO is always dissociated, 
whereas the longest period Cepheids have lower maximum temperatures and so never dissociate all 
their CO, thus have the largest color amplitude. 

The dependence of color on CO implies that the 4.5 /xm band should be subject to a metallicity 
effect, effectively ruling it out as a distance indicator. Fortunately, the bandhead is confined to 
the [4.5] filter alone; no CO or other molecular feature is present at 3.6 //m, and so it can still be 
confidently used to measure Cepheid distances. The physics behind the cyclical CO variation, along 
with the effect of metallicity on Cepheid's IRAC colors will be described in detail in Scowcroft et 
al. (2012, in preparation). 



6.3. PL Slopes and CO Absorption 

We now return to consideration of the slopes of the PL at different wavelengths . Figure [8] shows 



the de pendence of the slope on wavelength, a plot that is repeated from figure 4 of lFreedman et al 



(|2008l ). The [3.6] and [4.5] slope values have been updated using the Table 3 values in this paper. 
The slope of the [4.5] PL relation is clearly deviant. We can again understand this as due to CO 
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as follows: as we move from shorter to longer periods we are moving to cooler stars; these have 
more CO absorption on average. More CO means depressed 4.5 fim flux and thus a shallower PL 
slope. Therefore, this plot gives further evidence that the color curves are due to CO absorption 
that affects only [4.5]; the [3.6] magnitude is unaffected. 

A smaller magnitude effect is also seen in the K band where a shallower CO bandhead is 
present. However, this effect can be mitigated using the K$ band, which excludes the longest 
wavelength section of the bandpass — the part that happens to include the CO bandhead. 

7. Summary 

We have used the IRAC camera on "Warm Spitzer" to measure [3.6] and [4.5] light curves for 
85 Cepheids spread across the face of the LMC and covering a 6 to 140 day period range. Each 
star has been measured at 24 phase points. Our data agree well with the random phase data of the 
SAGE project, the latter taken for different purposes. 

The light curves (shown in Figure [2]) allow us to measure the mean magnitudes and color of each 
Cepheid to a precision never achieved before in the mid-infrared. These lead to our best PL relations 
based on the measurements of the [3.6] and [4.5] magnitudes for 67 stars with 10 < P < 60 days. 
The PL relations have been measured to such high precision that their scatter directly correlates 
with the front-to-back geometry of the sample of Cepheids in the LMC. This correlation recalls 
a similar effect presented in P04. Using the residuals of the period-luminosity relation, we have 
measured the tilt of the galaxy as 28°, with respect to the plane of the sky. 

The dense and uniform sampling of our Cepheid light curves reveals a cyclical color variation 
which we interpret as being due to the effect of temperature-sensitive molecular absorption of CO 
in the Cepheid's atmosphere. When the Cepheid is cool, the CO present in its atmosphere absorbs 
at 4.67 fim, suppressing flux exclusively in the IRAC 4.5 /jm band, making the star appear blue 
in the [3.6] — [4.5] color. As the Cepheid contracts, and so heats up, the CO molecules begin to 
dissociate, ending the suppression of the 4.5/iin flux and reddening the color. The temperature 
dependence of the CO absorption also explains the period-color relation we observe (see Figure EJ): 
the cooler, longer period Cepheids appear bluer (more highly absorbed by CO at 4.5 (j,m) than their 
short period counterparts. The implications of this temperature-sensitive effect will be studied in 
detail in Scowcroft et al. (2012, in preparation). 

The period-luminosity relations measured here will be used in combination with those from 
our Galactic sample (Monson et al. 2011, in preparation), to establish the bottom rung of the 
extragalactic distance ladder. They will then be combined with parallel observations of Cepheids in 
the Local Group and beyond, to establish the distance scale to a systematic accuracy not achievable 
before now. We will also be testing for the effect of metallicity on the mid-IR period-color relation, 
by combining this sample with Galactic and SMC Cepheids, (Scowcroft et al. 2012, in preparation). 
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9. Appendix 



The observing plan we adopted using Spitzer to optimally follow up known Cepheids is very 
differently constrained from what one would or could attempt from the ground. Ground-based 
observations are obviously restricted by 24-hour resonances with the day-night observing cycle. 
They are often additionally modulated by the lunar cycle (especially for optical observations). 
They are stochastically interrupted by weather, instrument availability, technical problems, and 
time assignments. Randomly spaced observations of known Cepheids may thus be the only way to 
fully sample light curves and to determine their mean values. In the case of space-based observations 
of known Cep heids randomly samp led data may be the only kind av ailable (e.g. data from the 
SAGE project (IMeixner et al.ll2006l ). which was mined for Cepheids bv lFreedman et al.l (|2008l ) and 
Madore et all l|200g| )). 



However, the 24-hour availability of a solar-orbiting telescope opens up a unique opportunity to 
design the observing program to be optimally scheduled on a star-by-star basis. Let us first consider 
the case of random sampling. In randomly sampling any given distribution function the error on the 
derived mean value goes down in inverse proportion to the square root of the number N of randomly 
timed observations (i.e., a oc 1/y/N). What is not widely appreciated is that convergence of the 
error on the mean can be very significantly accelerated by deterministic (non-random) sampling. In 
the case of an individual Cepheid where the time distribution of luminosities (i.e., its amplitude) is 
bounded (but not a priori known), and where th e cycling time (i.e., its peri od, but not necessarily 
its phase) is also known, then it is possible (see lMadore fc Freedmanll2005l ) to construct a simple 
observing strategy where the errors on the mean decrease in direct inverse proportion with the 
number of observations (i.e. a oc 1/N). Not surprisingly, for a specified number of observations N, 
and a known period of variation P, the optimal separation of data points in time is the uniform 
distribution P/(N+1). 

The distribution function in this case is closely approximated by a rectangular probability 
density distribution. The equivalent sigma for such a distribution is A/y/l2, where A is the total 
range of the distribution function; in the case of a Cepheid, A is its amplitude. In the mid-IR the 
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amplitudes of Cepheids can reach 0.8 mag (e.g. HV00883), which converts to an equivalent sigma 
of ±0.23 mag. If the required precision on a determination of the mean magnitude were to be set 
at 0.01 mag (i.e., one percent) then random sampling of the light curve would require over 500 
observations! To reach the same precision using deterministic (uniform) sampling a grand total of 
23 observations are required. To compensate for slight variations away from exact scheduling we 
requested and obtained 24 observations per LMC Cepheid. The light curves shown in Figure [2] are 
a testament to the care and efficiency with which Spitzer was scheduled for this program. 
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Fig. 1. — Photometry of two randomly chosen constant stars in the HV00872 (Star A) and HV02883 
(Star B) fields. The dashed lines show the average magnitude of the star and ±lcr limits. The 
dispersion in both cases is approximately 0.01 mag. The time baselines differ as the Cepheids in 
the fields were observed with different cadences. No systematic trends are seen. 
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Fig. 2. — IRAC light curves of the LMC Cepheid sample, in order of decreasing period. Point sizes 
are comparable to the uncertainties in the [3.6] (top) and [4.5] (middle) panels. The bottom panels 
show the variation of the IRAC color with phase. Upwards (more negative) in color corresponds 
to greater CO absorption as discussed in Section EJ 
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Fig. 3. — LMC mid-infrared PL relations. The solid lines indicate the relations given in Table El 
cut at 10 and 60 days, i.e. the open circles are not included in the fit. The dashed lines show ±2cr. 
Error bars are consistent with the size of the points. 
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Fig. 4. — Magnitude residuals from the [3.6] and [4.5] period-luminosity relations, with short period 
(P < 10 days) Cepheids as red open stars, long period (P > 100 days) as blue open triangles, and 
the PL fitting sample as black filled circles. The inset plot shows the whole sample; the main plot is 
cut at ±0.2 mag to show an expanded view of the details. The dashed line indicates A [3.6] = A [4.5], 
the solid line is a fit to the black points. The deviation from a one-to-one fit is believed to be due 
to CO pushing the [4.5] points further from the adopted PL. 
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Fig. 5. — Color and magnitude residuals from the period luminosity relation. The dashed lines show 
the la dispersion. Excluding the longest period (logP > 2) stars, the color dispersion around the 
[3.6] ridge line is 0.018 mag. Assuming that the instability strip has a regularly filled, rectangular 
distribution, the width will be y/l2a. This gives the color width of the instability strip as 0.062 
mag. 
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Fig. 6. — Positions of the Cepheids within the LMC. Open circles show Cepheids which lie above 
the PL relation (i.e. are brighter than average), filled circles show Cepheids which lie below it. 
The bottom left plot shows their positions relative to the center of the LMC, with the coordinate 
system rotated such that the line of nodes is horizontal. The top left and bottom right plots show 
the deviations from the PL relation separated into the £ and rj components. The top right plot 
shows the histogram of the size of the PL deviations. The dispersion in the PL is minimized by 
assuming a tilt of 28° in rj relative to the plane of the sky. 
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Fig. 7. — Period-color relation for our LMC sample. The solid line shows the PC relation of 
Equation[3j with the dashed lines showing Ha. The trend apparently reverses for the very longest 
period (logP > 2) objects. 



-32- 



B V R I 



J H K [3.6] [4.5] [5.8] [8.0] 



"i 1 1 1 — i — r 



"i 1 1 1 1 — i — r 



3.5 



Period-Radius relation 



a 

o 

GO 



2.5 



J I I I I L 



Madore & Freedman 1998 
Persson et al. 2004 
This Work 
Madore et al. 2009 

I I I I I I L 



10 



log A(ycxm) 



Fig. 8. — Relationship betwe en effective waveleng t h and the slope of the Cepheid PL relation. 
The PL sl opes are taken from Madore & Freedman! (|l998h (BVRI), P04 (JHK), this work ([3.6], 
[4.5]) and iMadore et al.l (|2009l ) ([5.8], [8.0]). There is an asymptotic trend towards —3.45, which is 
consistent with the theoretical prediction of a slope of —3.4 for the period-radius relation. However, 
the [4.5] (and to a smaller extent, K) slope is clearly discrepant. We attribute this to the CO 
bandhead that suppresses flux in the 4.5 fim. The [3.6] slope is not affected — it follows the trend 
of the other wavelengths — and is therefore appropriate to use for distance determinations. 
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Table 1. Magnitudes for LMC Cepheids. 
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Table 1 — Continued 
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14.466 


12.205 


0. 


.005 


12.221 


o. 


.004 


-0.017 


o. 


.002 


HV02337 


6.863 


13.247 


0. 


.016 


13.212 


o. 


.015 


0.034 


0. 


.005 


HV02338 


42.200 


in r on 

10.580 


0. 


nnT 

.007 


1 n f* A n 

10.640 


o. 


.007 


n a r ^7 

-0.057 


0. 


nno 

.002 


HV02339 


13.880 


12.139 


0. 


.004 


12.143 


o. 


.004 


-0.005 


0. 


.001 


HV02352 


13.632 


12.282 


0. 


.003 


12.267 





,003 


0.014 


o. 


.001 


HV02369 


48.377 


10.297 


0. 


.006 


10.355 


o. 


.006 


-0.058 


0. 


.003 


HV02405 


6.924 


13.309 


0. 


Oil 


13.283 


o. 


.012 


0.025 


0. 


.005 


HV02432 


10.918 


12.446 


0. 


.015 


12.427 


0. 


.016 


0.018 


0. 


.005 


HV02447 


118.350 


9.440 


0. 


.006 


9.438 


0. 


.005 


0.003 


0. 


,002 


HV02463 


13.965 


12.142 


0. 


.004 


12.134 


0. 


.004 


0.006 


0. 


,001 


HV02527 


12.950 


12.452 


0. 


.005 


12.456 


0. 


.005 


-0.004 


0. 


.002 


HV02538 


13.869 


12.246 


0. 


.003 


12.266 


0. 


.003 


-0.020 


0. 


.002 


HV02549 


16.220 


11.885 


0. 


.004 


11.891 


0. 


.004 


-0.006 


0. 


.001 


HV02579 


13.428 


12.121 


0. 


.003 


12.114 


0. 


.004 


0.006 


0. 


.001 


HV02580 


16.918 


11.879 


0. 


.006 


11.910 


0. 


.005 


-0.031 


0. 


.001 


HV02733 


8.722 


12.887 


0. 


.011 


12.855 


0. 


,010 


0.029 


0. 


.005 


HV02749 


23.105 


11.625 


0. 


.004 


11.682 


0. 


.004 


-0.058 


0. 


.002 
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Table 1 — Continued 



Cephcid 


Period 


[3.6] 


°"[3 61 a 


[4.5] 


°"[4 51 a 


[3.6] - [4.5] 


a \3 61 -14 5i a 




(days) 


(mag) 




(mag) 




(mag) 




HV02793 


19.220 


11.734 


0.006 


11.767 


0.005 


-0.033 


0.002 


HV02827 


78.802 


9.744 


0.004 


9.824 


0.003 


-0.081 


0.002 


HV02836 


17.528 


11.988 


0.005 


12.008 


0.005 


-0.019 


0.002 


HV02854 


8.635 


12.810 


0.016 


12.785 


0.013 


0.024 


0.005 


HV02883 


108.970 


9.953 


0.010 


9.958 


0.009 


-0.006 


0.002 


HV05497 


99.100 


9.345 


0.004 


9.376 


0.003 


-0.030 


0.002 


HV05655 


14.211 


12.308 


0.005 


12.325 


0.005 


-0.015 


0.002 


HV06065 


6.838 


13.301 


0.011 


13.282 


0.013 


0.014 


0.005 


HV06098 


24.237 


11.224 


0.003 


11.218 


0.003 


0.006 


0.001 


HV08036 


28.369 


11.329 


0.006 


11.374 


0.006 


-0.044 


0.002 


HV12452 


8.736 


12.837 


0.016 


12.832 


0.015 


0.003 


0.005 


HV12471 


15.863 


12.178 


0.004 


12.208 


0.004 


-0.028 


0.002 


HV12505 


14.389 


12.327 


0.005 


12.338 


0.005 


-0.013 


0.001 


HV12656 


13.400 


12.269 


0.002 


12.248 


0.002 


0.020 


0.001 


HV12700 


8.153 


12.975 


0.011 


12.953 


0.012 


0.022 


0.005 


HV12717 


8.843 


12.877 


0.015 


12.842 


0.015 


0.034 


0.005 


HV12724 


13.744 


12.419 


0.004 


12.444 


0.005 


-0.025 


0.002 


HV12815 


26.135 


11.212 


0.006 


11.257 


0.006 


-0.042 


0.002 


HV12816 


9.108 


12.890 


0.011 


12.869 


0.011 


0.019 


0.002 


HV13048 


6.853 


13.165 


0.011 


13.152 


0.013 


0.010 


0.005 


P71-U01 


22.558 


11.661 


0.007 


11.702 


0.007 


-0.038 


0.002 



a For Cepheids with P < 12 days the uncertainties scale with 1/y/N rather than 
1/N as the observations of short period Cepheids were not phase locked. 
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Table 2. Photometry of HV00872. 



MJD b 


[3.6] b 


a 


[3.6] 


[4.5] b 


o 


'[4.5] 






(mag) 






(mag) 






55133 


.665 


11. 


.538 


0. 


.013 


11 


.633 


0. 


.015 


55134 


.756 


11. 


.568 


0. 


.013 


11 


.605 


0. 


.014 


55136. 


.253 


11. 


.510 


0. 


.013 


11 


.498 


0. 


.013 


55136 


.953 


11. 


.485 


0. 


.012 


11 


.471 


0. 


.014 


55138 


.621 


11. 


.396 


0. 


.013 


11 


.372 


0. 


.013 


55139 


.561 


11. 


.378 


0. 


.012 


11 


.330 


0. 


.013 


55140. 


.963 


11. 


.334 


0. 


.011 


11 


.296 


0. 


.013 


55141 


.986 


11. 


.268 


0. 


.011 


11 


.276 


0. 


.012 


55143. 


.049 


11. 


.244 


0. 


.012 


11 


.255 


0. 


.012 


55144. 


.443 


11. 


.205 


0. 


Oil 


11 


.235 


0. 


.013 


55146 


.159 


11. 


.182 


0. 


Oil 


11 


.232 


o. 


.012 


55146. 


.916 


11. 


.162 


0. 


Oil 


11 


.235 


o. 


.012 


00140. 


.o4U 


11 


.165 


0. 


Oil 


11 


.243 


o. 


.012 


55149. 


.495 


11. 


.167 


0. 


Oil 


11 


.267 


o. 


,014 


55151. 


.036 


11. 


.171 


0. 


Oil 


11 


.290 


o. 


.012 


55151. 


.989 


11. 


.187 


0. 


Oil 


11 


.292 


o. 


.012 


55153 


.002 


11. 


.220 


0. 


Oil 


11 


.315 


o. 


.013 


55155. 


.321 


11. 


.266 


0. 


Oil 


11 


.394 


o. 


.013 


55155. 


.906 


11. 


.282 


0. 


Oil 


11 


.378 


o. 


.013 


55156. 


.963 


11. 


.302 


0. 


Oil 


11 


.419 


o. 


.013 


55158. 


.231 


11. 


.357 


0. 


.012 


11 


.460 


0. 


.013 


55159. 


.318 


11. 


.389 


0. 


.012 


11 


.528 


0. 


.014 


55160 


.793 


11. 


.439 


0. 


.012 


11 


.523 


0. 


.014 


55161. 


.870 


11. 


.472 


0. 


.012 


11 


.584 


0. 


.014 



Note. — The Cepheid photometry is avail- 
able in the electronic edition of the Astrophys- 
ical Journal. 

a MJD = JD - 2,400,000.5 d 

b Each magnitude point corresponds to a 
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flux weighted average of the five dither po- 
sitions. 
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Table 3: Period-luminosity relations 



Band 


Sample 


N 


a 


cr a 


b 




a 


[3.6] 


lOd < P < 60d 


67 


-3.41 


0.08 


12.74 


0.03 


0.108 




SAGE a 


70 


-3.40 


0.07 


12.73 


0.02 


0.135 




No cut 


85 


-3.16 


0.05 


12.67 


0.02 


0.134 




6 d < P < 60 d 


82 


-3.31 


0.05 


12.74 


0.01 


0.105 




P > 10 d 


72 


-3.14 


0.07 


12.67 


0.03 


0.142 


[4.5] 


lOd < P < 60d 


67 


-3.31 


0.08 


12.74 


0.03 


0.115 




SAGE a 


70 


-3.35 


0.07 


12.77 


0.03 


0.141 




No cut 


85 


-3.10 


0.05 


12.67 


0.02 


0.131 




6 d < P < 60 d 


82 


-3.22 


0.05 


12.70 


0.02 


0.112 




P > 10 d 


72 


-3.11 


0.07 


12.68 


0.03 


0.139 



a The SAGE data from 



Madore et al. 



(2009) were re-fit to be consistent with the CHP analysis. 



